Quantum Pumping with Ultracold Atoms on Microchips: Fermions versus Bosons 
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We present a design for simulating quantum pumping of electrons in a mesoscopic circuit with 
ultra-cold atoms in a micro-magnetic chip trap. We calculate theoretical results for quantum pump- 
ing of both bosons and fermions, identifying differences and common features, including geometric 
behavior and resonance transmission. We analyze the feasibility of experiments with bosonic 87 Rb 
and fermionic 40 K atoms with an emphasis on reliable atomic current measurements. 

PACS numbers: 37.10.Gh,03.75.Nt, 03.75.Ss,73.23.-b 



Pumping is any cyclical time-varying mechanism that 
generates sustained flow. Quantum pumping [l|, 0] in 
mesoscopic solid state circuits is a coherent quantum pro- 
cess for generating directed transport of charge with time 
dependent potentials, but no applied bias field. With 
its promise of precise and reversible flow control at the 
single electron level and extension to transport of spin 
3] and entangled electron pairs 0], quantum pumping 
has been the subject of considerable theoretical research 
Despite potential technological applications, quan- 
tum pumping experiments in solid state system have not 
been successful, partly due to dominant competing rec- 
tification effects associated with electrically charged car- 
riers 0, 0, Q • Neutral ultra-cold atomic systems present 
a possible path around the current bottleneck by avoid- 
ing such complications. An atomic circuit using a Bosc- 
Einstein condensate (BEC) or a degenerate Fermi gas 
(DFG) can test basic theoretical predictions, while also 
providing a reference for experiments in solid state sys- 
tems. 

In this letter, we present a design for an experiment to 
test quantum pumping theory with ultra-cold atoms in a 
micro-magnetic potential on a chip. Ultra-cold atoms 
open up the possibility of studying not only fermion 
quantum pumping but also boson pumping, as well as 
the influence of variable interactions and long range or- 
der, in a fully controlled and tunable system. We present 
theoretical results for both types of atoms in prototypical 
pumping schemes and we analyze the feasibility of a cold 
atom based experiment with numerical simulations. 

Mesoscopic circuit with atom chips: A prototypi- 
cal mesoscopic circuit consists of a device, e.g. a quantum 
dot, connected by nanowires to macroscopic contacts. At 
low temperatures, electrons and holes can have mean free 
paths longer than the nanowires, so they can be described 
as freely propagating particles in one dimension (ID). 
The device presents a scattering potential for the parti- 
cles, so that transport is reduced to a scattering problem 

We can simulate this setup with the atom chip based 
scheme shown in Fig. [TJ Atom chips are substrates on 
which currents in lithographically imprinted wires gen- 



erate a micro-magnetic trapping potential for ultra-cold 
atoms. These chips can efficiently produce both BECs 
and DFGs with temperatures in the 10 nK to 1 /iK range 

[m eh. Two reservoirs connected by a ID quantum 

wire can be implemented on an atom chip by two 3D 
micro-magnetic traps connected by a quasi-lD magnetic 
guide, generated by co-propagating currents in two paral- 
lel wires (red in Fig.[T]) on the substrate, with a constric- 
tion for the tighter ID section. The atoms are trapped in 
the plane of the wires, with the substrate between them 
removed [l2[ , which also allows optical access from above 
and below. The trapping potential is harmonic along all 
principal axes, including the axial one due to a current 
through the two 'end cap' wires (dashed-green in Fig. 1) 
below the trapping plane [131 ] . Residual defects in the 
trap potential can be suppressed by applying an AC cur- 
rent through the principal trapping wires, while keeping 
the external axial magnetic field and the current in the 
end-cap wires constant 1J] . The "device" , or scattering 



potential, can be realized with a dipole laser focused onto 
the ID section. 




FIG. 1: (a) Configuration for generating two micro-magnetic 
trap reservoirs connected by a 1-d channel. The red wires 
provide radial confinement, while the dashed green 'end-cap' 
wires, located 50 /im below, provide axial confinement. The 
large red arrows are probe lasers for measurements on the 
trapped atoms, represented by the blue structure. The ver- 
tical (purple) laser implements the pump potential; (b) a 1 
fiK equipotential for alkali atoms trapped by 250 mA and 10 
mA in the red and green wires, respectively, along with a 1 
G axial magnetic field (the transverse/axial scale is 37); (c) 
transverse isopotential curves along the 1-d channel from 50 
to 1000 fiK showing its symmetry and significant trap-depth. 
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The generated atomic current can be determined from 
a measurement of the momentum distribution of the par- 
ticle flow, since the average current can be written as 
J = (h/m) J \ip(k)\ 2 kdk. Bragg spectroscopy [l5[ is ide- 
ally suited for measurements of the momentum distribu- 
tion, since it can be selectively applied to atoms in the 
1-d channel and combined with fluorescence imaging for 
high signal-to-noise detection. A spectroscopic flag can 
be attached to the kicked atoms by adding the hypcrfinc 
splitting to the base detuning of the probe lasers, thus 
changing their hyperfine level. A large fraction of the 
scattered photons can be collected by a microscope lens 
located a few millimeters above the atom chip and im- 
aged onto a high sensitivity camera. We calculate that 
roughly a hundred photons per atom can be detected 
with a fluorescence pulse of a few hundred microseconds. 

Theory of bosonic and fermionic pumps: Quan- 
tum pumping has been studied exclusively for fermions in 
solid state systems, and primarily in the adiabatic regime 
where the pump period exceeds the dwell time of the car- 
riers at the potential. With atomic experiments in mind, 
we extend the theory of pumping to include bosons. 

As with electrons in nanowires, the dynamics of 
atoms in the central segment is quasi-lD with quasi- 
continuum description along the transport direction 
and quantized transverse channels (n). The ax- 
ial and transverse components can be factorized [lij ]. 
^{x,y,z,6;i) = £ n ip n (z,t)<f> n (r), (in cylindrical sym- 
metry); J dz\ip n (z, t)\ is the population fraction in the 
n-th channel. Scattering influences the evolution of the 
axial functions, with little effect on the transverse pro- 
files. For weak interactions 171 ] . phase fluctuations of 
degenerate bosons can be neglected in the ID section, 
so the axial dynamics has an effective description in 
terms of a ID non-linear Schrodinger equation (NLSE) 
5 [-d* + (z/7) 2 ] ^ + giD\ip\ 2 ip = id t ip where 7 = uj r /u; z 
is the aspect ratio. Fermions and non-interacting bosons 
are both described by setting g\ E = 0. The axial po- 
tential variation is small over the ID segment, so we set 
z 2 /j 2 ~ 0, allowing a plane wave description. The radial 
trap frequency uj r sets the scale for our expressions: the 
energy, length and time units are hoj r , l r =-\/fi/(niuj r ) and 
lu^ 1 . Time evolution is governed by the axial energy of 
available channels E=/i — nTioj r , where fi is the chemical 
potential . 

For a slowly varying external potential V(x, t), an adi- 
abatic expansion to first order in the time-derivative ap- 
proximates the time-dependent scattering states 

V> fc ,„0M) ^ <»(a:) - iJdx'G t (x,x';E)d t 7Pl n (x'), (1) 

in terms of the solutions of the Lippmann-Schwinger (LS) 
equation |^ n ) = \<f>k,n) + G t {E)V t \4>k,n) for the instan- 
taneous potential V t {x). RereG^E) = [E-Hq-V*]- 1 is 
the instantaneous Green's operator and \4>k,n) are scat- 
tering states of the time- independent Hamiltonian Hq, 



and k,n label wavevector and channel. The zeroeth or- 
der current vanishes for a pump with no bias field jij. 
Denoting the second term in Eq. (Q]) by Atfj k ^ rn (x, t) , the 
adiabatic pumped current of spin polarized fermions is 

J F (x,t)=j: n JdEf(E)J^6(^-E)x (2) 

Im{ip k * n (x)d x A^ n {x,t) + A^ n (x,t)d x ^l n (x)} 
= 27 jrdEf(E)d E Jdx'V(x , ,t)lm{G t ^(x',x)d x G t E (x,x')} 

using Green's function identities and the LS equation 0]. 
At low temperatures, the Fermi distribution is approxi- 
mately a step function, f(E) ~ 6(Ep — E), so we obtain 



Jp(x, t) 



dx' ■ 

— V(x')lm{G* EF (x', x)d x G Ep (a;, x')} (3) 

Z7T 



In order to describe bosons, we use the parabolic 
dispersion of plane waves to replace J^°dEf(E)-^ — > 

/ °° dkf(^-)jr Tz >° j Q dk-^j: . For non- interacting degen- 
erate bosons at rest or in a superposition of momentum 
states, I A;) and | — fc), the pumped current is 

J B (x,t) = ^Jdx'V(x')lm{G E (x',x)d x G E (x 1 x')} (4) 

Thus the pumped current for a degenerate fermi gas at 
fermi- vector k is related to the pumped current for degen- 
erate bosons at wavevector k by Js(fc) = 27r^ Jp (fc 2 /2). 
For a BEC at rest, one sets k = 0, after the derivative. 

Essential features of quantum pumps can be under- 
stood with models involving time-varying single barrier 
potentials: with variable strength Vi(x,t) — U{t)5{x) or 
translating uniformly V2 (x, t) = US(x—vt). For adiabatic 
variation, Eqs. ([3]) and (T?|) give the pumped currents: 



J { p\x,t) = sgn(x) T 



Ukj 



(2) _ 1 k F vU 



(U 2 + k 2 



(U 2 + k 2 ) 2 



They show the role of symmetry: V\ generates no net 
particle transport from one reservoir to the other over a 
period, due to antisymmetry with respect to x; while V% 
being symmetric leads to net transport. In general, both 
symmetric and antisymmetric parts can be present. 

The fermionic current for V2 is always in the direction 
of motion of the potential, but the bosonic current can 
flow opposite ( Fig. [2]): When the bosons have sufficient 
energy k 2 > U, the transmitted fraction dominates, and 
particles going against the barrier have a higher transmis- 
sion probability; for fermions, the averaging over states 
washes out this effect. Over a period T, the net pumped 
particles, J x T, is independent of the velocity v and de- 
pends only on the parameter path traversed by the poten- 
tial. It is a geometric quantity analogous to a geometric 
phase (l8j , a feature shared by all adiabatic quantum 
pumps. 
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Bosons 



Fermions 



FIG. 2: Pumped current versus velocity of a uniformly trans- 
lating potential: (a,b) <5-barrier, (straight lines are adiabatic 
approximations); (c,d) square barrier of width a = 10 l r and 
height U/a. Barrier strengths and wavevectors, (U,k) for 
bosons and (t/, ftp) for fermions, are shown. 

The pumped current at arbitrary barrier velocities for 
Vi can be found using a Galilean transformation to be: 



Je 



2vU 2 {U 2 +v 2 



k 2 ) 



U 4 +2U 2 (v 2 +k 2 ) + (v 2 ~ k 2 ) 2 



Jf- 



■In 



U 2 +(k F +v) 2 



U 2 +(k f 



(6) 



The adiabatic expressions in Eq. ([5]) are retrieved with 
a Taylor expansion for v -C k, hp. The fermionic current 
vanishes for kF=0, as the number density vanishes; but 
for a stationary BEC, Js(fc=0) = vU 2 /(v 2 + U 2 ), en- 
tirely due to reflection. At high barrier velocity v k, 
Jb ~ J B (k=0) and J F ~ k F vU 2 /(i:(v 2 + U 2 )) like the 
adiabatic limit with kp and v interchanged. 

Even for a, finite translating square barrier (SB) V$b — 
(U /a)8(x — vt)6(—x + vt + a) of width a and height U/a, 
analytical expressions for the pumped current can be like- 
wise calculated, too lengthy to be shown, but plots based 
upon those solutions are shown in Fig. [2j The pumped 
current differs dramatically from the case of the delta 
barrier: (i) the finite height allows particle energy to ex- 
ceed the barrier potential leading to sharp transitions at 
i(fc±t>) 2 = U, the classical cutoffs for transmission, and 
(ii) the finite width creates oscillations due to resonance 
transmission. For bosons (Fig. [2fc)), the oscillations are 
pronounced, with the current vanishing and reversing for 
some velocities; but less so for fermions (Fig. [2](cL)) due 
to averaging over wavevectors. For a translating barrier, 
both classical and quantum features are manifest, but for 
a translating well, the behavior is quantum mechanical. 

The quantum nature of pumping becomes truly signif- 
icant in a turnstile pump comprised of two barriers with 
heights oscillating out of phase with each other. This 
model has been studied [g, \l9^ for fermions, and here 
we present results for bosons contrasted with fermions. 
Essential features can be understood with two delta func- 
tion potentials U±(t)5(x =p a) with oscillating strengths 
U + (t)=l + cos{tut) and U- (t)=l + sin(wi), that trace out 
a circle over a period T=2tt/lo. In this limiting case, 
the current is entirely due to quantum interference [l9| . 
Reversing the cycle reverses the flow. 

The insets in Figure [3] show that the currents on the 
left and the right of the potential are not in sync and vary 
over time, but their time integrals over a full cycle are 




FIG. 3: A double barrier pump of bosons (a)-(c) and fermions 
(d)-(f). Particle transport in a pump cycle as a function of: 
(a,d) barrier separation, 2a, at k—kp—1; (b,e) wavevector k 
(fermi-vector kF for fermions) at o=l; (c,f) angular frequency 
u>, at a=k=kF=l- Insets show current versus time at the right 
and left reservoirs . The length and time units are l r and ■ 

equal. Resonant transmission effects are prominent due 
to finite barrier separation, 2a: The particle transport, q , 
in a pump cycle displays oscillations and peaks as a func- 
tion of the barrier separation, Fig. [3ja,d) , and also as a 
function of the wavevector k or kp, Fig.[3Jb,e). Fermions 
display less pronounced resonance behavior, due to av- 
eraging over momentum states. There has been recent 
interests in testin g re sonance transport through double 
barrier structures [201 ] , quantum pumps demonstrate this 
by periodic cycling of the potentials. The geometric na- 
ture of adiabatic pumps is clearly seen in Fig.[3fc,f), since 
particle transport per cycle is independent of uj. 

Feasibility Analysis: Pump potentials can be im- 
plemented with blue-detuned lasers at 532 nm focused 
to 1-5 /im gaussian-profile barriers. The lasers need to 
translate at velocities v ~ l r uj r ~ 0.5 cm/s or vary in 
intensity at frequencies ~ ui r ~ 2 kHz, easily achievable. 

Bosonic pumps at non-zero |fc| can be implemented 
with a broad (relative to pump potential) wavepacket 
split into counterpropagating momentum states by a 
Bragg pulse [Hj]. In this scenario, the reservoirs can be 
removed. For 87 Rb in the F — 2, raj — +2 state in 
the set-up of Fig. [TJ the transverse and axial trap fre- 
quencies are u r .iD = 2n x 5.1 kHz and u> a xiai — 2tt x 3.6 
Hz. For a wavepacket of 100 87 Rb atoms with scattering 
length of a s = 99gq a variational calculation 2l| yields 
the effective ID non-linear constant g\o = 67.3 and ax- 
ial Thomas-Fermi width 587 l r . After the axial trap is 
turned off and the Bragg pulse applied, the split wave 
packet evolves in the presence of the pump potential. 
We simulate this numerically by solving the NLSE with 
a split-step operator method; results are shown in Fig. 4. 
In the absence of nonlinearity (gid — 0) , the wavepacket 
simulations with square barriers are consistent with an- 
alytical results (Figs. 2(c) and 3(b)) obtained assuming 
plane waves, validating the method. Figure [4] also shows 
that Gaussian profile barriers and nonlinearity lead to 
some qualitative changes, but the pumped current or 
charge remains significant. The nonlinearity reduces the 
signal somewhat, and for the turnstile, the broader bar- 
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FIG. 4: Numerical simulation of pumping with interacting 
and non-interacting bosons. In (a) k=l. Square barrier (SB) 
plots show consistency with analytical results. The inset com- 
pares simulations for a SB turnstile of separation 2a=2 with 
analytical results for a 5— barrier turnstile as in Fig. 3(b). 
As the barrier width decreases, with area fixed at 1, the SB 
results approach the 8~ barrier analytical results. 

riers and barrier-separation lead to more closely spaced 
oscillations in the current as packet velocity (k) varies. 
Numerical simulations [221 ] show that the pump signal 
for the turnstile is more sensitive to chip trap roughness 
than the translating barrier scheme, but AC suppression 
of roughness 14[ is sufficient for a robust signal. 

A 1000 atom wavepacket has initial peak density 
7.3 x 10 14 cm~ 3 and chemical potential fi^d = 0.26/j.K = 
1.09?kiV- The number of atoms can be significantly in- 
creased, considering: (i) More atoms mean stronger non- 
linearity and faster expansion, requiring longer traps to 
allow sufficient interaction times with the pump; with- 
out the reservoirs the axial length can be extended upto 
1000 fim. (ii) To remain in the transverse ground state 
(for single channel), /i < 2hu r ; our variational calcula- 
tions gives /i ~ 1.6fud r with Nin = 2.0 x 10 4 atoms. 

For fermion pumps with 40 K the currents listed for 
Fig. [1] produce trap frequencies a factor of yrn^JrriK — 
1.5 higher than with 87 Rb. Energetically, the ID section 
can contain 0J r .iD/^>axiai — 1400 spin-polarized fermions 
in the lowest transverse channel due to the Pauli prin- 
ciple. Since the size of harmonic oscillator eigenstates 
scales as V2N, for the axial oscillator length l z = 6.9 /im, 
the 250 /xm ID section will hold about 700 atoms; each 
reservoir contains 50 times more. The lowest channel can 
accommodate Fermi vectors up to kp = 1.4Z" 1 . 

Conclusions: Our analysis has shown that quan- 
tum pumping experiments can be done with current 
atom-chip technology, allowing a broad survey of a 
process that has eluded confirmation in solid state 
systems. In addition to simulating fermion pump- 
ing, ultra-cold atom based experiments open up the 
possibility of studying quantum pumping of bosons 
which we expect to show enhanced resonant tunneling 
and current reversal. Furthermore, the scheme can 
be adapted to search for conductance quantization 



for pumping with periodic lattices [l[ by imposing a 
moving optical lattice on the 1-d quantum channel. In 
a broader context, our design is easily adapted to a 
variety of mesocopic transport experiments, important 
in electronic systems, like conductance quantization and 
spin transport, yet hardly explored with ultracold atoms. 
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